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I . NOMENCLATURE 



Cross-sectional area of inlet to fuel grain 
Cross-sectional area of fuel grain 
Mass transfer parameter, (h^ - h Q ) /h Q - h^) 
Coefficients in k-e turbulence model 
Specific heat at constant pressure 
Friction coefficient without mass addition 
Inside diameter of fuel grain 
Generation term in equations for k and e , 




Mass flux of air based on bulk velocity 
Enthalpy, C p (T - T°) + m Qx H Q /i 
Effective heat of combustion, cal/gm of fuel 
Stoichiometric ratio of oxygen mass to fuel mass 
Turbulence kinetic energy 

"Length scale" of eddies at near wall grid point P 
Length of fuel grain 
Mass fraction 

Mass flux from fuel surface 
Average molecular weight 
Pressure 

Molecular Prandtl number 
Radial coordinate 
Regression rate of fuel surface 
Universal gas constant 

Near wall grid Reynolds number, p V y r) /l J 0 
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*p 


= Internal radius of fuel port 


R t 


1/2 

= Turbulent Reynolds number, y p P/U Q 


s * 


= Source tern in differential equation for generalized dependent 
variable <p 


T 


= Temperature 


V 


= Velocity 


y P 

+ 

y P 


= Distance from wall to near wall grid point P 

1/2 

= Non-dimensionalized distance from wall, (p/y ) (x /p) y„ 

O o r 


z 


= Axial coordinate 


Y 


= Slip function from Couette flow for boundary condition on k 


£ 


= Rate of dissipation of turbulence energy 


n 


= Combustion Efficiency based on total temperature rise across 
combustor 


V 


= Molecular viscosity 


U eff 


= Turbulent viscosity 


p 


- Density 


°eff 


= Effective turbulent Prandtl or Schmidt number 


T 


= Shear stress 


<p 


= Generalized dependent variable 


ip 


= Stream function 


(1) : 


= Vorticity 



Subscripts 



fu 


= Fuel 


h 


= Enthalpy 


n 2 


= Nitrogen 


o,w = 


= Wall 


of : 


= Refers to the "property" m, - m /i 

fu ox 
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ox 



= Oxygen 

P = Near wall grid point 

r = Radial component 

T = Reference condition deep within fuel grain 

2 = Axial component 

0 - Tangential component 

Superscripts 

o = Reference conditions 

- Average value 
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II. INTRODUCTION 



Several characteristics of the solid fuel ramjet indicate that it may be 
superior to other forms of propulsion for tactical weapons used at interme- 
diate ranges and high speed. Having no moving parts, the solid fuel ramjet 
is simple and relatively inexpensive to fabricate. Weight of the system and 
its threat as a fire hazard are both decreased by use of the solid fuel. 

To be used in a tactical situation, the solid fuel ramjet has to demon- 
strate combustion stability and efficiency over the expected operating 
envelope of altitudes and Mach numbers. It must also show performance 
comparable to that of liquid fuel ramjets and ducted rockets. 

Combustion studies on the solid fuel ramjet have been underway at United 
Technologies-Chemical Systems Division since 1971. Initial work showed low 
temperature rise combustion efficiencies. The discovery of the rearward 
facing step at the combustor entrance as a flameholder was a significant gain 
in solid fuel ramjet technology. Overall performance, however, was reduced 
by the high stagnation pressure loss produced by the rearward facing step. 
Further work by United Technologies in the field of flame stabilization in- 
volved various inlet designs including aerogrids, distorted flows, non- 
circular inlets and vortex generators. This effort led to the development of 
inlets which minimized the required inlet step height and decreased the 
effects of inlet distortion. 

Work was also underway at the Naval Postgraduate School on the internal 
ballistics of the solid fuel ramjet. Early work‘d showed that inlet turbulence 
and distortion may have a significant effect on flame stability. 

The solid fuel ramjet, which uses air as the oxidizing agent, is similar 
to the hybrid rocket and has two distinct combustion zones within the fuel 
grain (Fig. 1) . Behind the step is the recirculation zone where an intense 
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mixing of reactants and products takes place. The hot products ignite the 
reactant and combustion may in the limiting case approach that of a well- 
stirred reactor. This combustion region acts as the flame initiator for the 
combustion which takes place further down the fuel grain. 

Downstream of the flow reattachment point a boundary layer develops. 

Here combustion is similar to that of the hybrid rocket. A diffusion flame 
exists in the developing turbulent boundary layer between the fuel rich zone 
near the wall and the oxygen rich central air core. Heat is transported by 
convection and radiation to the solid surface which causes decomposition of 
the fuel. Studies have shown that the rate of decomposition of the fuel is 
sensitive to the combustion pressure, the inlet air temperature and the mass 
flux of the air. Finite rate kinetics may be of importance for certain 
fuels and operating environments. 

The combustion process in the solid fuel ramjet is thought to be mixing 
limited at combustion pressures greater than 10-15 psia for all-hydrocarbon 
fuels. However, this remains to be verified and other fuels have demon- 
strated behavior more characteristic of kinetically controlled combustion. 

Unbumed gaseous fuel escapes from under the flame at the aft end of the 
fuel grain and results in decreased combustion efficiency. This resulted in 
the use of a mixing chamber downstream of the fuel grain. Temperature rise 
efficiencies have been found to be a function of the rate of mixing between 
the fuel rich boundary layer and the cental air core. Too rapid mixing may 
quench the chemical reactions and further reduce combustion efficiency. Two 
somewhat distinct combustion zones also exist in the aft mixing section, the 
recirculation zone formed by the fuel grain or aft orifice and the central 
core zone G 
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The use of aft-end mixing devices to improve combustion efficiency is an 
area of recent interest in solid fuel ramjet technology. As another means 
of promoting mixing, the bypassing of a portion of the inlet air around the 
fuel grain and dumping it into the aft mixing chamber is being re-evaluated. 

A meteorlogical sounding rocket capable of an altitude of 200,000 feet was 

2 

designed and built by Anderson, Greenwood and Company, Houston, Texas in 1961 . 
The MET JET employed a ramjet using a magnesium and magnesium-aluminum alloy 
epoxy-metal charge. Eighty-five percent of the inlet air was bypassed around 
the fuel grain to mix in an afterburner section with the fuel rich primary 
flow. 

In the bypass systems the flow rates into the fuel grain inlet and aft 
mixing section are critical factors in determining combustion efficiency. 

Bypass air flow of too low a percentage of the total air mass flow or of too 
low momentum may have negligible effect on the combustion efficiency. Bypass 
air of too high a flow rate or momentum may have a negative effect on the 
combustion process. The bypass air is of appreciably lower temperature than 
the species which exit from the fuel grain, and air injected at too high a 
rate may cool the process sufficiently to affect the kinetics of the reaction. 
Furthermore, the combustion process in the aft mixing chamber may also be a 
function of the axial and radial positions and the angular orientations of 
the aft dumps. 

This experimental investigation considered the effects of bypass con- 
figuration and operating conditions on the combustion behavior of a solid fuel 
ramjet which utilized polymethylmethacrylate as the fuel. 

In order to guide development efforts and to provide needed input to 
systems analysis studies, combustion models are needed which can predict the 
fuel regression rate, flammability limits and combustion efficiency as a 
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function of hardware design and operating environment. 

Modeling efforts for the combustion process in solid fuel ramjets haye 

3 

proceeded along two directions. Dunlap has developed a semi-empirical 

regression rate model which is based upon reacting turbulent boundary layer 

theory and experimentally obtained wall heat transfer rates in non-reacting 

1 4 

flows. Work at the Naval Postgraduate School * has been directed toward 
adaptation of the Spalding, et al model'** ^ for heat and mass transfer in 
recirculating flows to the SFRJ geometry. 

Experimental data ^ have shown that the recirculation zone is fuel rich 
and that the primary combustion region (or flame zone) spreads out from along 
the shear layer between the recirculation zone and the inlet flow. The 
flame has been observed to be quite broad near the flow reattachment posi- 
tion and subsequently develops into a turbulent diffusion flame within the 
developing boundary layer downstream of reattachment. 

The initial modeling work at the Naval Postgraduate School considered 
the effects of finite rate kinetics and combustor geometry on the flame 
pattern and internal flow field. 

Several weaknesses were evident in the earlier work. Fuel regression 
rates were not calculated. Also, the flame pattern was not in agreement 
with experimental data, i.e., the recirculation region was predicted to be 
oxidizer rich and the flame within the boundary layer region spread too 
rapidly toward the centerline of the motor. In addition, the model had 
not been adequately checked for accuracy against experimental data for 
velocity, pressure and turbulence intensity. 

The present investigation was conducted to improve the qualitative 
accuracy of the model, to incorporate prediction of the fuel regression 
rate, and to check the validity of the model against experimental data. 
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The majority of the analytical studies were conducted for polymethyl- 
methacrylate (PMM) fuel. PMM-air combustion produces minimal radiative heat 
transfer to the fuel grain for high mass flow rates through the port. The 
latter condition is characteristic of non-bypass systems. The model 
developed for this system did not include radiative heat transfer or finite 
rate kinetics. 

The current investigation also considered the application of the model 
to the combustion behavior of all-hydrocarbon fuels which exhibit significant 
gas-phase radiation. In particular, model predictions were compared to the 
data obtained by Schadow \ 

IIIo EXPERIMENTAL INVESTIGATION 
A. Equipment and Procedures 

A schematic of the solid fuel ramjet is shown in Fig. 2. The apparatus 
employed a fuel port to dump inlet area ratio of 9.0 and a fuel port to 
nozzle throat area ratio of 4.0. The exit area of the grain (for most tests 
conducted) was held fixed at the initial port area by using a thin orifice 
plate. The aft mixing chamber had a length to diameter ratio (L/D) of 2.93. 
The aft mixing chamber consisted of three interchangeable sections such that 
the axial location of the bypass dump could be varied. One section contained 
four inlets for the bypass air. The inlet diameters were 2.70 cm with plugs 
available to reduce the diameters to 2.065, 1.908, 1.458, or 0.635 cm, or to 
seal the inlets off entirely. An additional section was fabricated which pro- 
vided bypass dumps oriented tangentially to the mixing chamber wall and 
perpendicular to its centerline. In most tests conducted the bypass dump was 
located in the forward position, within the recirculation zone of the aft 
mixing chamber. 
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All test firings were performed in the jet engine test cell at the Naval 
Postgraduate School, Data for calculating temperature rise efficiencies based 
on combustion pressure were obtained while varying primary and bypass air flow 
rates and bypass dump geometry. Nominal air flow rates were 0,1 and 0.2 lbm/ 
sec for non-bypass runs and 0.2 lbm/sec for runs using bypass. A summary of 
the test conditions are presented in Table I. 

Table I Summary of Experimental Conditions 

Air Inlet Temperature (°K) : 292, 389 

Bypass Dump Location: forward and aft of flow reattachment 

Bypass (% Port/% Bypass): 100/0, 65/35, 50/50. 35/65 

Momentum Ratio (Dump /Port) 

for 50/50 Bypass: 0.3, 0.5, 0.6, 1.0, 5.3 

Dump Configuration: none, 2 at 180°, 2 at 90°, 4 at 90°, 2 at 180° with swirl 

The majority of the test firings were made using a nominal inlet air 
temperature of 70°F. For those using the air heater (non-vitiated air), two 
to three hours were allowed for the temperature to stabilize at the ramjet 
inlet and in the bypass air line. 

The weight, length and inside diameter of the fuel grains were measured 
prior to being mounted in the motor. The ignition sequence normally lasted 
for six seconds. An average of three seconds was required for the ignition 
flame to propagate from the head-end assembly into the fuel grain. After 
five seconds, the air was directed through the ramjet motor. Ignition was 
continued into the first second of the run to insure that combustion would 
be sustained. Total time of the ignition flame in contact with the fuel 
grain amounted to approximately three seconds. Two tests were made using 
only the ignition system in order to determine the rate of consumption of 
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the PMM grains during the oxygen-methane ignition. These data were used to 
correct the initial weight of the fuel used in the efficiency calculation. 

Combustion normally lasted for forty-five seconds. The motor was ex- 
tinguished at the end of each run by simultaneously venting the air to the 
atmosphere and actuating the nitrogen purge system. Low pressure air was 
then blown through the motor for cooling. 

Thirty-three hot firing PMM tests were conducted. Combustion efficien- 
cies were based on the total temperature rise from just upstream of the rear- 
ward facing step at the fuel port inlet to the nozzle inlet. Inlet total 
temperature was derived from the measured inlet static temperature. Actual 
combustion temperature was derived from the measured static pressure in the 
aft mixer 0 

Regression rate was calculated based on weight loss of the fuel. The 
inside diameter of the aft end of the fuel grain was also measured before and 

o 

after each firing. However, as observed by Boaz and Netzer °, it was found 
that weight loss gave a more consistent value of regression rate than the 
method based on aft end diameter change, due to the non-uniform regression 
along the length of the grain. 

The calculated temperature rise efficiencies based on pressure are prone 
to error, since uncertainty in efficiency varies as the square of the uncer- 
tainty in the pressure measurement. Errors as large as 7% were possible. 
However, the uncertainty in percent theoretical C* was approximately 2% and 
the temperature rise efficiencies followed the C* efficiencies within 1.5%. 
B. Results and Discussion 

(a) Regression Rate 

Several runs were made without using bypass air. These showed the depen- 
dence of the regression rate, r (in/sec) on combustion pressure, P (psia) 
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and average mass flux of air, 



G . 
air 



(lbm/in -sec) 



to be: 



r = 0.0043P ,29 G . .38 
air 



( 1 ) 



A plot of the regression rate versus this empirical regression rate equation 
for the 0.2 and 0.1 lbm/sec non-bypass runs made with cold inlet air is 
shown in Fig. 3. 

g 

In their work with PMM, Boaz and Netzer had shown a regression rate 
dependence of the form: 



= CP’ 51 T 



.34, 



.41 



air 



air 



( 2 ) 



Insufficient runs were made in this investigation at other than ambient 
inlet temperature to calculate a temperature dependence for the regression 
rate. The above equations for regression rate agree closely for the depen- 
dence on G. The dependence on pressure was significantly less in this study. 

In the work of Boaz and Netzer, chamber pressure was varied intention- 
ally from 37 to 108 psia by varying the throat size. In this study pressure 
varied between 33 and 63 psia only as a result of varying combustion effi- 
ciency and G. The dependence of regression rate on G results from convective 
heat transfer to the fuel surface whereas the dependence on the pressure can 
result from radiation and/or finite rate kinetics. 

With the application of bypass air, the dependence of regression rate on 
pressure and air flux was altered. In Fig. 4 are plotted the 0.2 lbm/sec 
non-bypass cases and the bypass runs made with 0.2 lbm/sec total air flow rate 
at the nominal test condition. The nominal test condition was a bypass air 
flow rate of 0.1 lbm/sec and two dumps with diameters of 2.065 cm. This 
corresponded to a dump momentum to fuel grain port momentum ratio of approxi- 
mately 0.5. A slightly stronger dependence on pressure was shown while 
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regression rate indicated very little or no dependence on air flux for the 
values tested. Here regression rate took the form: 

r = 0.00116P* 42 G . •° 03 (3) 

In the bypass situation, the mass flux through the grain is low but the 
pressure is maintained high due to the total mass flux through the nozzle 
throat. However, correcting the regression rate for the increased pressure 
per equation (2) does not result in regression rates as high as the experi- 
mental data. In addition, regression rates based on weight and diameter 
agreed, indicating that the change did not result from different combustion 
behavior within the aft mixing chamber. These conditions of low G and high 
P minimize the convective heat flux. At lower G the regression rate 
increases relative to the air flux, i.e., more unburned fuel exists under- 
neath the diffusion flame within the thicker boundary layer; thus more gas 
with radiative properties is present. 

(b) Combustion Efficiency 

In Fig. 5 are plotted the cases from both Figs. 3 and 4. As can be 
seen from this figure the use of bypass air flow in a solid fuel ramjet using 
polymethylmethacrylate as a fuel has the effect of reducing combustion 
efficiency. For the non-bypass runs, a decrease in air flow also causes a 
slight decrease in combustion efficiency. While maintaining the same air 
flux through the grain, injecting bypass air into the mixing section brings 
a further decrease in performance. Decreasing air flux through the fuel 
grain while maintaining the same total flow rate in the bypass case also 
brings about a reduction in combustion efficiency. Decreasing the air flux 
through the fuel grain increases the percentage of fuel that must be burned 
in the aft mixing chamber. Apparently, the mixing chamber was not of 
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sufficient length and/or the bypass air caused quenching of the mixing 
chamber combustion. 

To further study the decrease in performance found when using bypass, 
other test conditions were considered. Two runs were made with the dumps 
located behind the aft mixing chamber reattachment point. There was only a 
slight decrease in combustion efficiency (81%) over that of identical runs 
in the forward position (83%). 

Three runs were made varying the momentum of the individual dumps from 
the nominal test condition. Small changes in momentum had no effect. Vary- 
ing the momentum ratio from 0.3 to 1.0 did not significantly change the effi- 
ciency. However, when a significantly large increase in momentum was 
effected (5.3), the combustion performance showed a noticeable decrease (73%). 

Only limited data were taken for other than two bypass dumps at 180 
degrees. However, the data showed that with low dump momentum two or four 
dumps spaced at 90 degrees reduced the combustion efficiency. Four dumps at 
ninety degrees with high momentum decreased the combustion efficiency slightly. 

The last variation in geometry studied was the use of swirl in the bypass 
dump process. This was accomplished by injecting the bypass air with a tan- 
gential velocity component. The solid fuel ramjet motor was also run without 
fuel grain ignition, both with and without swirl, to determine the effect of 
the induced vorticity on the effective exit nozzle diameter. Decreasing the 
effective nozzle diameter would increase the combustion chamber pressure and 
give false indications of higher efficiency. The swirl was found to not 
affect the effective throat diameter. Two test firing runs were made with 
swirl, giving an average combustion efficiency of 86%, the highest perfor- 
mance of any bypass run. 
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Every form of bypass used caused a decrease in combustion efficiency. 

This indicates that the decrease in performance is due to the effect of the 
bypass air flow on the kinetics of the combustion process within the aft 
mixing chamber. Essentially equal decreases in performance were obtained 
with the dump air injected both behind and in front of the reattachment point. 
Apparently a significant portion of the total combustion process takes place 
downstream of the reattachment point in the aft mixing chamber. This agress 
with the temperature data presented by Schadow ^ for an all-hydrocarbon fuel. 

In the case of PMM, the light weight unburned hydrocarbons that enter the aft 
mixing chamber apparently burn most completely when allowed to react slowly 
with the available oxygen in the hot flow from the core of the fuel grain. 

The possibility of the temperature of the air entering the aft mixer 
causing reduced combustion efficiency was considered. A run was made in which 
the temperature entering both the grain and aft mixer was increased from 
approximately 292°K (520°R) to 389°K (750°R). The combustion efficiency did 
not change significantly, indicating that for PMM-air combustion the dump air 
temperature is probably not as important as the quantity and temperature of 
the unburned fuel and the mixing rate. 

In the region in front of the reattachment point, the less that was done 
to disturb the flow resulted in better performance. In the case of swirl 
where the bypass flow remained close to the wall, bypass efficiency was maxi- 
mized. When the bypass momentum was increased to where the flow disturbed the 
fuel rich layer between the recirculation region and the air rich central core, 
performance decreased. This again indicates that a major portion of the com- 
bustion process takes place along this fuel rich layer or that the process 
downstream is highly dependent on the high temperatures in this layer. It 
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also indicates that the combustion mechanisms around this layer are more 
important than those which occur within the recirculation region. These re- 
sults imply that different optimum bypass dump configurations and momentum 
should be expected for different fuel systems. 

It is known that the use of bypass improves performance in solid fuel 
ramjets using all-hydrocarbon fuels. In the case of oxygen-containing fuels 
and for fuels which decompose into monomers or small hydrocarbon molecules 
(such as PMM) , results indicate that the fuel burns most efficiently without 
bypass . 

The use of bypass systems has meant an increase in weight, cost and com- 
plexity of the solid fuel ramjet. In addition, they may introduce combustor- 
feed system coupling. The use of a fuel which has sufficient density impulse, 
regression rate and flammability limits to minimize inlet total pressure 
losses has led to the use of all-hydrocarbon fuels. Although PMM does not 
meet the criteria for a good fuel, the results of this study indicate that 
future fuel studies may be fruitful if directed toward ones which contain low 
percentages of oxidizer and/or substances which unzip the hydrocarbon chain. 

(c) Combustion Pressure Oscillations 

In the non-bypass test runs, the inlet and combustion pressures exhibited 
a steady, small amplitude (approximately 2% of chamber pressure) oscillation 
of approximately 150 Hz. In the bypass runs at nominal conditions, the same 
oscillation appeared though of considerably higher amplitude (approximately 
30% of chamber pressure). 

In the case of low or very high aft dump momentum, a second oscillation 
appeared along with that previously mentioned. This oscillation was of very 
low frequency (1 Hz.) and large amplitude (approximately 20% of chamber pres- 
sure) and may be connected with behavior within the mixing chamber recirculation 
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zone. The bypass-generated very low frequency pressure oscillations are in 
the range of frequencies which could be coupled to the fuel regression rate. 
Other fuel/feed systems may exhibit similar or different oscillatory behaviors 
depending upon the feed system design and the kinetics of the combustion 
process in the aft mixing chamber. 

C. Conclusions and Current Work 

The use of bypass air flow in a solid fuel ramjet which utilizes PMM as 
the fuel decreases the combustion efficiency. Air flow injected into the aft 
mixing chamber has a more pronounced effect on the combustion process when a 
high enough momentum is provided for the bypass air to reach the fuel rich 
shear layer trailing from the port of the fuel grain. 

A significant amount of combustion occurs downstream of the reattachment 
point in the aft mixing chamber. Bypass air injected into the aft mixing 
chamber also has an effect on the regression rate upstream in the fuel grain. 
For PMM fuel burned at high pressures and low air mass flux, the principal 
mechanism for wall heat flux became radiation, and resulted in the regression 
rate becoming insensitive to the air flux through the fuel port. 

With the use of bypass systems it is possible to set up combustor-feed 
system type oscillations and instabilities dependent upon the effects of the 
bypass flow on the aft mixing chamber combustion process and on the fuel 
regression rate. 

Current experimental work is being directed at the effects of bypass on 
the combustion behavior of all-hydrocarbon fuels. In addition, the effects 
of design and test conditions on combustion pressure oscillations are being 
investigated. 
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IV. ANALYTICAL INVESTIGATION 



A. Model Overview 

The basic assumptions were that the flow was steady, recirculating, 
two-dimensional and subsonic. In addition kinetic heating and vorticity 
sources resulting from spatial gradients in the effective viscosity were 
neglected. The turbulent Lewis number was taken to be unity. The trans- 
formed variables of vorticity (eo/r) and stream function (¥) were used 
together with enthalpy (h) , and species mass fractions (oxidizer, fuel, 
etc.). A modified Jones-Launder^ * ® two equation turbulence model was 
used to calculate the viscosity throughout the flow field. Thus, two addi- 
tional variables must be considered, the turbulence kinetic energy (k) and 
the turbulence energy dissipation rate (e) . The modified Jones-Launder 
model incorporates five empirical constants and the effective viscosity is 
calculated using the following expression: 



In the present study the combustion was considered to be mixing 
limited with a one-step simple chemical reaction. Thus, 1 gm fuel + i 
gm oxidizer (1+i) gm of products and fuel and oxidizer could not exist 

together. Species considered were oxygen, nitrogen, fuel and products. 

The present form of the model neglected radiation and density was calculated 
using p = pM/RT . 

With these assumptions the governing equations for axisymmetric flow 
can be cast in the standard elliptic format^* 




( 4 ) 




convection 

terms 



diffusion 

terms 



source 

term 
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where the parameters are presented in Table II. Additional discussion of 
these equations and utilization of various turbulence models can be found 
in Ref. 11 and 12. 



Table II. 



Equation Parameters 



5,6 




a 



h 



a r = a„ - a. = 1.0, a =1.3 
of N 0 k e 
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B . Boundary Conditions And Solution Procedure 

The inlet flow was considered to be "plug flow” although other inlet 
conditions (distortion, etc,) are easily incorporated. The inlet turbu- 
lence intensity was taken to be 5,2% in order to closely match experimental 
measurements. Axial gradients for all variables were considered zero at 

the flow exit and were taken as zero for all "conserved properties” 
m 

( m fu ~ , m^ , h) at the step-face. Radial gradients for all varia- 

bles except ¥ were taken as zero at the axis of symmetry. The major 
area of concern in the specification of boundary conditions for the SFRJ is 
along the fuel surface. When calculations are made with the model for an 
adiabatic or isothermal wall with no mass addition, inaccuracies in the 
turbulence model are somewhat maskedo The fuel mass addition at the wall 
and the subsequent combustion greatly amplify the need for accurate effec- 
tive viscosity calculations near the wall if observed experimental behavior 
is to be approximated. 

All "conserved” properties when in dimensionless form, had identical 
governing equations and boundary conditions. Thus, only one equation was re- 
quired to be solved for all "conserved” properties. The resulting set of five 

0 ) 

equations (— , ip, k , £, m^ - m^/i) were cast in finite difference form and 
solved using the Gauss-Seidel method with upwind differences and relaxation. 
Under relaxation was required: 0.5 for ^ and ip , 0.3 for k and e , and 
0 o 9 for m^ - m /i • A non-unif ormily spaced 17 x 25 rectangular grid system 
was employed. Solution was performed on an IBM 360-67 using Fortran H and 
required approximately thirty-five minutes to complete 1000 iterations. 

The fuel surface was considered to have a uniform and constant temperature 
and the mass fraction of oxidizer at the surface was considered to be zero. 

In order to obtain a fuel rich recirculation zone, a boundary layer flame 
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pattern in agreement with experiment, and accurate conservation of species, 
several modifications had to be made to the original model^*^. For simpli- 
city the boundary layer was considered to consist of two parts, a laminar 
sublayer and a turbulent layer. The border between the two regions was 
considered to occur where y^ = 24 . At each near wall grid point y"*" was 
calculated. 



If y + > 24: C, /2 « 0.02846/R °* 25 

p fo e 



<«) 



13 

The velocity profile was corrected for blowing and y __ was calculated 

ef f 

using Eq. (4). 



If y < 24: C. /2 = 2/R . 

J V — fo e 



(7) 



The laminar viscosity was calculated using 

y = y (T/T ) 0,5 
o o 



( 8 ) 



In the actual calculation procedure, the grid was selected to insure 
that y^ _< 24 . The value of 0.02846 is required for matching of the 
two- zone boundary layer. 

The boundary condition for vorticity was calculated from 
T 2 

(-) = — - — where x = (C, /2) p? In (1 + B) /B (9) 

r o u r o fo 

o w 

in which the wall shear stress is corrected for blowing by using the simple 
Couette flow approximation ln(l + B) /B . The control volumes used at the 
fuel surface for the micro-integration are shown in Fig. 6. The control 
volume for all variables but k and £ were taken to the wall. For k 
and £ the control volumes were taken to the half-grid point. Gradients 
at the wall for "conserved" properties were taken to be of the form 
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(10) 



I 



(|~) rt = {(h - h )/y } In (1 + B) /B 

dy o p op 



The wall mass flux provided the boundary condition for ¥ and was 
calculated from 



m” = (— ) A /(h - h ) = (— ) ln (1+3) 

p 3y o v o T 'p y 

r r r ■’p 



The boundary condition on k was calculated from 



( 11 ) 



k Q = Y k p (12) 

Y = - 1.0 , R < 19.8 

= - 0.39, R > 19.8 
t 



The boundary condition on £ was considered to occur at the near 
wall grid point 



where 



- 1C 3/2 M d 


(13) 


p p 


= 0.2 y , y + < 24 
P P 


(14) 


= 0.4 y , y + > 24 
P P 


(15) 



The boundary conditions for — , k and £ on the step face were 

calculated in a similar manner except that there was no wall mass addition, 

none of the micro- integration control volumes were extended to the wall, and 

l = 2.44 y 
P P 

14 

Launder and Spalding have discussed the strengths and weaknesses of 
the k-£ model using both the wall function and low-Reynolds-number modeling 
methods. They indicate that modification of the boundary conditions on £ 
were sometimes found necessary on the step face in order to obtain results in 
agreement with experiment. In this investigation wall-functions were employed 
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for the k-e model. Since the mass addition rates were low (B < 1.5) the 

five constants in the model^ , ^ > ^^' were not modified. When y + < 24 along 

P 

the fuel grain the effective viscosities predicted using the k-£ model were 
not used at the near wall grid points. In this case, the boundary conditions 
given by Eqns. (13), (14) and (15) were required to provide realistic values of 
^eff an< * fraction at the second grid point from the wall. An alter- 

nate procedure could have been to use the low-Reynolds-number modeling 
methods^* ^ In recent work with smaller inlet step heights it has been found 
necessary to employ Eqns. (13), (14), and (15) for £^ along the step face. 
Much additional work is apparently required to better define the near-wall 
region. 

The finite-difference equations^’ ^ require the use of space averaged 
values for the source terms presented in Table II. At the near wall grid 
points large gradients exist if y^ + > 24 . As mentioned above, this was 
prevented at the fuel surface by appropriate selection of the radial grid 
spacing. Along the step face, however, integrated average values had to be 
calculated for — and all velocity gradients. Velocity profiles along the 

step face for y + > 24 were assumed to obey the 1/7 th power law. 

P 

In order to obtain realistic solutions the stream function at the grid 
point downstream of the step inlet was required to be held fixed at the 
value which existed at the step corner. 

C. Results and Discussion 

Divergence would occur in the iteration unless the initial conditions 
specified for oo/r and ip were in close agreement (for example, 1/7 th 
power law distributions). It was also necessary to prescribe an initial 
length scale distribution. The length scale distribution was used to cal- 
culate initial values of £ from k using Eqn (13) . The £ field was 
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% 



held constant for an initial 35 iterations in order to promote convergence. 
Specification of the initial length scale distribution required trial and 
and error methods for other than the simple rearward facing step geometry. 
Variations in gas density were also prevented from varying for the first 50 
iterations. 

At each axial location species and energy conservation calculations 
were made. False diffusion errors were larger in the recirculation zone 
where the grid lines were at large angles to the streamlines. Small errors 
in the solution of the "conserved" property equation presented some difficul- 



ties since 



m > , m « m and m_ _ 

fu 9 ox 9 pr Nr 



were all calculated from the one 



pr » 2 

solution. Calculation of fuel conservation was especially difficult since 



fuel concentration gradients near the wall were very large. In general, 
energy and species (other than fuel) were conserved to within five percent. 

With these accuracies, fractional changes between iterations were a maximum 
of 0.3 percent. This required between 750 and 1000 iterations. The "conserved* 



OJ 



properties and stream function converged more rapidly than did , k and e , 



CjO 



For the higher inlet velocities, the relaxation parameter for — had to be 



reduced to values between 0.3 and 0.4. 

The location of the "reattachment point* 1 was found to be somewhat sen- 
sitive to the distance from the step face to the first axial grid line 
This resulted from specification of the stream function as unchanging as 
discussed above. However, once a particular grid was chosen which yielded 
a reattachment point in agreement with experiment for one step height and 
inlet velocity, good results were obtained for other step heights and inlet 
conditions. The reattachment point was found to be insensitive to inlet 
velocity and was found to move slightly upstream with low rates of wall mass 



addition. Both of these results are in good agreement with experimental 
8,15 



data 



26 



A true reattachment point does not exist when mass addition through 
the wall is present. However, a pseudo-reattachment point was calculated 
by assuming it to occur where the average velocity was zero between the wall 
and the near-wall node. A comparison of this ff reattachment point” with 
chemical reaction is made with experimental data from non-reacting flows in 
Fig. 7. It should be noted that the predicted flow reattachment occurred 
downstream of the Krall-Sparrow^ maximum heat transfer point but upstream 
of the Phaneuf-Netzer^ data with 10% wall mass addition. 

Typical predicted streamline patterns and temperature distributions 
are presented in Figure 8 for plexiglas fuel and air. The dividing stream- 
line is also sketched in the temperature field. The recirculation region 
was calculated to be fuel rich and the flame initiated along the shear layer 
and spread out as it approached the reattachment region. The flame region 
was located within the developing turbulent boundary layer downstream of 
the reattachment zone. These observations are in qualitative agreement with 
experiment. 

A comparison of predicted plexiglas fuel regression rates with experi- 
ment is presented in Figure 9. It should be noted that the location of the 
maximum regression rate obtained experimentally was affected by the inlet 
air distortion/ turbulence intensity. The model predicts a maximum regression 
rate upstream of the experimental values and slightly upstream of the Krall- 



age fuel regression rate was in good agreement with experiment and the regres- 
sion pattern was qualitatively correct (although shifted upstream) . In the 
present model only convective transfer to the fuel surface is considered. The 
predicted regression profile may not be correct for fuel systems which pro- 
duce large amounts of radiative heat transfer to the surface. 




location for maximum heat transfer in non-reacting flow. The aver- 
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Figure 10 compares the predicted centerline turbulence intensity 
(assuming isotropic turbulence) with experimental data for non-reacting flow 
The behavior with axial distance was generally in good agreement but the 
inlet region exhibited behavior not observed experimentally. The model 
overestimates the velocity increase of the air as it enters the combustor, 
which in turn reduces the predicted turbulence intensity. 

Figure 11 compares experimental and predicted axial pressure profiles. 

The experiment and theory were not for identical conditions, but in general 
it is seen that the model predicts the pressure to level-off upstream from 
experimental measurements. This was expected since the model underestimated 
the distance to the reattachment point* Pressure calculations were very 
sensitive to solution accuracy since they depend upon second derivatives of 
stream-function. 

Figure 12 presents predicted radial temperature profiles as a function of 
inlet step size and air flow rate at a location downstream of flow reattach- 
ment. As the air flow rate was reduced for a fixed step size, the peak 
temperature moved further from the wall, the core flow became hotter, and 
the overall fuel-air ratio increased. This results from the fact that regres- 



sion rate varies as the 0.3 to 0.5 power of G a ^ r • Actually the model predicted 



that i ~ G^_!^ » whereas experimentally it was found for plexiglas that r ~ 



The recirculation zone in the solid fuel ramjet acts as the flame sta- 
bilization device. It operates fuel rich but, when operating within certain 
limits, provides enough energy to heat the air in the shear layer to the 
point where the downstream boundary layer flame can be stabilized in 
the high velocity flow. As the inlet flow rate was increased the tempera- 
ture in the upstream portion of the recirculation zone was found to increase. 
This resulted from the increased air flux into the fuel rich portion of the 



28 



zone. However, this same increased air flux into the recirculation zone 
decreased the temperature in the downstream portion of the zone and at the 
reattachment position. Thus, although the peak temperatures in the recircula- 
tion zone increased with increased inlet air velocity, the oxygen reaching the 
downstream flame zone decreased in temperature. In addition, the boundary 
layer flame narrows and moves closer to the wall. Eventually a condition would 
be reached where the oxygen reaching the boundary layer flame region would be 
too low in temperature to sustain the narrowed flame which exists in the rela- 
tively high velocity region close to the fuel surface. Earlier work with the 
1 4 

model * , in which one-step finite rate kinetics where employed, also demon- 
strated that increased inlet velocity can cause flame blow-off in either the 
recirculation region or the boundary layer region. 

Decreasing the inlet step size decreased the length of the shear layer 
(reattachment occurred further upstream) and decreased the air entrainment 
into the recirculation region. Although the shear layer was shorter in 
length, contact time with the recirculation zone did not necessarily decrease 
because the inlet velocity was lower. However, the model predictions showed 
a decreased temperature in the recirculation zone, resulting from a still more 
fuel rich mixture, and lower temperatures near reattachment. Again it would 
be expected that a flammability limit would be reached with increasing inlet 
diameter. 

Figure 13 presents the calculated percentage of unburned fuel as a function 
of air flow rate. As discussed above, the fuel mass fraction was the least 
accurate calculation. Thus, the results can only be interpreted 
qualitative. As G ^ is decreased, an increasing amount of unreacted and 
hotter fuel is passed into the aft mixing region downstream from the end of 
the fuel grain. 
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Schadow has made measurements of temperature and species concentrations 
in a reacting solid fuel ramjet . His experiments utilized an all-hydrocarbon 
fuel which results in significant radiative heat transfer to the fuel surface 
and higher fuel regression rates 0 However, the data are all that is currently 
available for SFRJ model validation in the reacting flow environment* Because 
the regression rate predicted using only convective heat transfer was too low 
for the radiative system, the predicted regression rates along the length of 
the fuel grain were arbitrarily increased by 30%. This cannot be expected to 
yield the correct axial variations in regression rate but was done to obtain 
a more realistic value of the blowing parameter for the higher regression 
rate system. 

Figs. 14, 15, and 16 present the theoretical temperature profiles compared 
to the data of Schadow^ for different air mass fluxes and axial locations 
within the fuel grain. The theoretical temperature profiles are not for 
exactly the same Lp/Dp locations as the experimental data. The profiles at the 
nearest corresponding grid line have been presented. 

Fig. 14 presents data for the high mass flux condition. The model pre- 
dicts more complete mixing and therefore more uniform temperatures in the 
recirculation zone (Lp/D„ = 0.52) than the data. In the region near reattach- 
ment (Lp/Dp = 2.3) the experimental data indicates an apparent ,! two-f lame 11 zone 
It is possible that just upstream of reattachment one flame may exist very 
close to the wall and another could result, farther from the wall, from the 
recirculating, reacting shear layer. However, the data of Schadow has not 
been corrected for radiation effects and an Lp/Dp of approximately 2.0 should 
be the most downstream location of reattachment for the inlet step height 
employed. If this "two-f lame’ 1 region is characteristic of the flow near 
reattachment, then the model fails to predict this detail. There could be 
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several reasons for this possible weakness in the model. Radiation is 
neglected in a system which produces large amounts of radiation. In addition 
the K-e model and/or the boundary conditions specified on the fuel surface 
for shear stress could result in a "smearing-out" of the species distribu- 
tions in the radial direction near the wall where large property gradients 
exist. Downstream of reattachment where the boundary layer develops 
(Lp/Dp = 5.2, 6.1) the profiles are in better agreement with experiment. The 
centerline temperature increase with axial distance is reasonably predicted. 

Similar results for lower air fluxes can be observed in Figs. 15 and 16. 
The model, without radiation-convection coupling to determine the fuel 
regression rate, predicted a stronger dependence of regression rate on air 
flux than obtained experimentally. This will also affect the predicted 
temperature profiles. The model does correctly predict the increasing 
boundary layer thickness with resulting peak temperatures further from the 
wall and higher centerline temperatures as the air flux is decreased. The 
regression rate variation with axial distance is different for the all- 
hydrocarbon system than for the convection dominated PMM fuel system. The 
model does a reasonable job of predicting this profile for PMM. For both 
fuels the regression rate increases from the head-end to near the reattachment 
zone and then decreases. Further downstream the PMM regression rates (and 
the model predictions) continue to decrease and/or level off whereas the 
all-hydrocarbon regression rates begin to increase again. 

The model assumed mixing limited and complete combustion so that 
"products" consisted only of H^O and CC^. This is obviously an oversimplifi- 
cation but perhaps warranted when compared to the other simplifications 
employed in the model v Schadow^ has measured with a probe the radial 
variation in moles of free 0^ per mole of at the exit plane of the fuel 
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grain. A comparison of his data with the model is presented in Fig. 17. 
Schadow^ found that in order to predict the correct variation in combustion 
temperature with mixture ratio in the fuel port, 40% of the carbon present 
had to be assumed as unreactive soot particles. As a qualitative attempt to 
compare the model (with infinite rate kinetics and no free carbon) with the 
data of Schadow, 40% of the carbon in "CC^ products” was assumed unreactive. 

The oxygen from the CO^ was combined with the unreacted 0^ to produce a new 
profile for the moles of ”free” 0^. The results are presented in Fig. 17. 
and show reasonable agreement with experiment. 

D. Conclusions and Current Work 

Temperature, pressure, turbulence intensity, and species profiles as well 
as the flow reattachment position were in qualitative agreement with experi- 
ment. The largest discrepancies existed in the recirculation region and 
near the fuel grain at reattachment. Regression rate profiles were in good 
agreement except for the all-hydrocarbon fuel far downstream from reattachment. 
The model appears to be adequate for qualitative evaluation of the effects of 
geometry and test conditions on the flame pattern and the amount of unburned 
fuel entering the aft mixing chamber. In order to obtain better quantitative 
results, radiation and finite rate, multiple-step kinetics should be incor- 
porated into the model. Further work is required to validate the k-£ model, 
especially in the near-wall regions. 

The HJ-0) variables make it difficult to obtain accurate predictions of 
the pressure field and to specify boundary conditions on the walls. The 
point iterative method employed also requires a large storage space and 
long computation times. 
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Current work with the model includes (a) extending its application to 



include the aft-mixing geometry as shown in Fig. 1 and (b) using the 
primary variables of pressure and velocity in place of the ¥-0) variables. 
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Fig. 1. Schematic of Solid Fuel Ramjet 
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Fig. 2. Schematic of Solid Fuel Ramjet Apparatus 
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Fig. 3. PMM Regression Rate- Non- Bypass 




Fig. 4. PMM Regression Rate, Bypass 
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Fig. 5. PMM Combustion Efficiency vs. Air Flux 
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Fig. 6. Near-Wall Control Volumes 
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Fig. 7. Reattachment Locations for Axisymmetric Flows 
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Fig. 8. Predicted Streamline and Temperature Distributions 
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Fig. 9. Plexiglas Regression Rates 




Fig. 10. Centerline Turbulence Intensity 
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Fig. 11. Axial Pressure Distributions 




Fig. 12. Predicted Radial Temperature Distributions 
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UNBURNED FUEL, 




Fig. 13. Predicted Fraction of Unburned Fuel at Grain Exit 
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Fig. 14. Comparison of Theory with Temperature Data of Schadow (G . = 34.5 

2 2 air 
gm/cm -sec (0.49 lbm/in -sec), P = 8.47 atm, D = 7.11 cm, 

L = 50.8 cm, T . = 605°K, A /A C = 2.74) P 

p air pi 






Fig. 15.. Comparison of Theory with Temperature Data of Schadow (G . = 16.2 

2 2 air 
gm/cm -sec (0.23 lbm/in -sec), P c = 6.73 atm, = 7.11 cm, 

L = 35.6 cm, T . = 590°K, A /A. = 3.55) 

p ’ air ’pi 
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Fig. 16. Comparison of Theory with Temperature Data of Schadow (G . =13.4 
2 2 air 

gm/cm - sec (0.19 lbm/in - sec), P = 6.8 atm, D = 7.11 cm, 

c p 

L = 50.8 cm, T . = 550°K, A /A. = 3.20) 

p ’ air ’pi 




Fig. 17. Comparison of Theory with Concentration Data of Schadow 
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